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. Several attempts to construct the thermodynamical formulation for non-equilibrium steady states 

' have been made along the context by Oono and Paniconi, called the steady state thermodynamics 

CN ; (SST). In this paper we study the SST for radiation field. Consider a cavity in a large black 

body, and put an another black body in it. Let the temperature of the outer black body be 
different from that of the inner one, and keep their temperatures unchanged. Then the radiation 
field between the black bodies are retained in a non-equilibrium steady state with a stationary 
energy flow. (The inner black body may be considered as a representative of a black hole with 
the Hawking radiation, and the cavity be the universe including the black hole.) According to the 
total momentum carried by photons par a unit time in passing through a virtual unit area, the 
pressure for the SST is defined. Then, taking this pressure as the basis, we propose the definitions 
■ of steady state thermodynamic functions and state variables satisfying the plausible properties: the 

O ' convexity /concavity of thermodynamic functions, the intensivity/extensivity of state variables, the 

Gibbs-Duhem relation and the Legendre transformation among thermodynamic functions. 

PACS numbers: 04.70.-s, 04.70.Dy, 05.70.-a, 05.70.Ln 
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I. INTRODUCTION 



along the context by Oono and Paniconi called the steady state thermodynamics (SST) 0. In this paper we 
O [ study the SST for radiation field. The reason we consider the radiation field is the future expected application to 
black hole physics. The black hole is the spacetime region of extremely strong gravity so that no matter can escape 
^ from the region in the context of classical physics (the theory of general relativity) ^3] • However it is well known that 
ly-j the quantum effect of the matter field on the black hole spacetime makes the black hole emit a thermal radiation 
00 from its event horizon (Hawking radiation) 0]. This means that the black hole can evaporate due to the energy loss 
by the Hawking radiation. While it is valid to assume that the black hole itself is in an equilibrium state, but the 
radiation field emitted outside the black hole should be considered in a non-equilibrium state with the energy flow 
by the Hawking radiation. However, so far, the matter fields surrounding the black hole have been assumed to be in 
, an equilibrium state, and it has been impossible to trace some details of the evaporation process with changing the 
I black hole's temperature and entropy at once :4j. Therefore we can not avoid considering a non-equilibrium state of 
^ ' the radiation field for a total understanding of the black hole evaporation. 

• Although our final aim is the application to black hole physics, we set the black hole physics aside and concentrate 
C on the system which can be designed for a laboratory experiment. Consider a cavity in a large black body, and 
^ , put an another black body in the cavity. (The inner one is a representative of the black hole with the Hawking 
^ ' radiation.) Let the temperature of the outer black body Tout be different from that of the inner one Tin{^ Tout), 
and keep their temperatures unchanged. Then the radiation field between the black bodies are retained in a non- 
equilibrium steady state with a stationary energy flow. Hereafter we abbreviate the non-equilibrium steady state 
5^ , with a stationary energy flow to the steady state. We treat this radiation field throughout this paper and search the 
consistent definitions of steady state thermodynamic functions and state variables for the radiation field within the 
phenomenological framework. 

Further we do not consider any statistical and quantum fluctuations in this paper, and assume that both of black 
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FIG. 1: Non-equilibrium steady state with a stationary energy flow for a radiation field. A black body of temperature Tin is 
put in the cavity made in a large black body of temperature Tout- The house- keeping heat rate Jhk retains the radiation field 
in a steady state. 

bodies emit always thermal radiations of their temperatures respectively without respect to the temperature difference. 
Then the energy Jhk which we should supply to the radiation field par a unit time in order to retain the steady state 
is given by 



where a is the Stefan-Boltzmann constant and Ain is the surface area of the inner black body. This Jhk is the 
house-keeping heat rate explained below. 

Along the context of the original consideration of the SST 0, we need to separate the total heat supplied to the 
radiation field during a change of steady state into two parts, house-keeping heat and excess heat. The house-keeping 
heat is the necessary heat supplied to the radiation field in order to make it stay at one steady state (without any 
change of the state). The longer the time interval for keeping the steady state becomes, the larger the house-keeping 
heat we should supply. So it is useful to consider the house-keeping heat rate suppHed par a unit time, which is 
given by eq. 1^. In changing the state of the radiation field from one steady state to another one, the total heat 
rate Jtot supplied to the radiation field at each moment during this change depends not only on the path of the 
change of steady state in the state space (of steady states), but also on the speed and acceleration of the change along 
the path in the state space, the geometrical shape of the system and the temperatures of black bodies. Therefore 
it is necessary to supply more (less) energy than Jhk to raise (lower) the temperature difference Tin — Tout, and 
consequently Jtot is generally different from Jhk at each moment during the change of steady state. Then the excess 
heat rate J^x is defined by Jex = Jtot — Jhk, which is the part of Jtot required to drive the change of steady state 
of the radiation field. Note that we find J ex = Jtot for equilibrium cases {Jhk = 0). The central spirit of the SST 
is the hypothesis that the excess heat in the SST plays the role of the total heat in equilibrium thermodynamics 
to define the SST entropy which specifies the reversible/irreversible processes for the steady states. Consequently 
we postulate that the phenomenological formulation of the steady states (SST) for the radiation field exists after 
extracting the effects of Jhk from various quantities, and that the SST has the same mathematical form as the four 
laws of the equilibrium thermodynamics using the excess heat in place of the total heat of equilibrium cases. Note 
that this postulation denotes the following three requirements. The first requirement is the existence of the state 
variables which characterise each steady state completely and uniquely. The second one is that the SST formulation 
includes the state variables describing the non-equilibrium order which vanish in equilibrium cases. Finally, the third 
requirement is to keep the basic properties of state variables unchanged: the convexity/concavity of thermodynamic 
functions, the intensivity/extensivity of state variables, the Gibbs-Duhem relation and the Legendre transformation 
among thermodynamic functions. 

In section m the definition of SST pressure are proposed, whose value is operationally determinable by laboratory 
experiments. Referring to this SST pressure in section lTTll we suggest the consistent definitions of SST free energy and 
the other state variables after introducing the intensive state variable specifying the non-equilibrium order. Section 
II VI is devoted to show the consistency of our definitions, but the SST second law is left to be confirmed in the future 
work. Summary and discussions are located in section Ivl 
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FIG. 2: The photon which contributes to Ap°"'. 
II. SST PRESSURE 

For constructing the SST for the radiation field, we should make state variables be determinable operationally. 
Therefore we search the definition of the pressure for the SST (SST pressure) using the total momentum carried by 
photons passing through a unit area par a unit time. Then an operational method for measuring the SST pressure 
is given according to the definition. But before considering the steady state, we start with reviewing the equilibrium 
case. 

A. For equilibrium case 

Let the temperatures of black bodies be equal Tin — T'ont = ^egi then the radiation field between them is in an 
equilibrium state of the pressure Peq{Teq) = (4fT/3c)Tg^, where c is the speed of light. In the region filled with the 
radiation field (between the black bodies) , image a virtual closed smooth surface C surrounding the inner black body. 
Then, as is well known, the relation Peq — Apj^ holds, where Apj^ is the total normal component of the momentum 
to C carried by photons passing through C par a unit area and a unit time. We divide it as Ap_L — Ap| + Ap|, 
where Ap| is the part due to the photons crossing C from the exterior region of C, and Ap| is due to the photons 
crossing C from the interior region of C. Note that, while all photons contributing to Ap| are emitted at the surface 
of the cavity, but the photons contributing to Ap| consist of ones emitted at the inner black body's surface and the 
others emitted at the cavity's surface. So further wc can divide as Ap| = Ap™ + Ap°"*, where Ap™ is the component 
emitted at the inner black body's surface and Ap""* is that emitted at the cavity's surface. Therefore we find 

Pe,-Apx + (Apf + Apf"*) . (2) 

Next let us place a very thin wall on C whose inner and outer surfaces are covered with perfect mirrors, and keep 
the temperature of the mirror wall be absolutely zero so that the mirror wall emits no photon and just reflects the 
photons emitted at the other black bodies. Then the radiation field in the exterior region of the mirror wall is in an 
equilibrium state with the pressure Pgq, and the radiation field in the interior region is also in the same state. For 
this case, as is well known, we find P^q = AJ| — A/|, where A/j^ is the normal component (to C) of the impact given 
to a unit area on the mirror wall par a unit time by the photons in the exterior region, and A/| is the impact by 
the photons in the interior region. Further, due to the perfectly elastic collision of photons at the mirror wall, the 
relations A/j^ — 2Ap| and A/| = 2Ap| hold. Consequently we obtain 

Ap; = Apf + Apf"*-iPe, (3) 

Here, after removing the mirror wall, let us calculate Ap°"* explicitly using the black body spectrum of photons 
emitted at the cavity's surface. Consider a photon emitted at a point q on the cavity's surface to a point r on the 
closed surface C as shown in figure |2] The magnitude of the momentum of this photon is hv/c, where uj is the 



4 



frequency of the photon. The normal component of this momentum to the surface C at r is (huj/c) cos 6, where 9 
is the angle between the momentum and the normal direction to C. Note that, because of the inner black body's 
shielding, no photon comes to r from the cavity's surface for some range of the angles 9 and (j), where 4> is the zenithal 
angle with respect to the normal direction to C at r. We introduce a function e{9, (p) where e = for the range of 
the angles 9 and with no coming photon to r and e = 1 for the other range of the angles. Using the black body 
spectrum, the number of photons coming to r from q par a unit area on C, a unit time, a unit frequency and a unit 
solid angle is given by 



cos 9 , 
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where sm9 arises with the integral about the solid angle around r, a; = jShw and a is the factor determined by 
the foUowings: the geometrical shape of the cavity and the inner black body, the configuration of the inner black 
body in the cavity, the geometrical shape of C and the choice of the point r on C From eq. (jSJ, we find 
Ap™ = Peql'^ — Ap°"*. This means that Ap™ has the complicated dependence on the geometry of the system as 
Api"* does. But note that eq. |(2Jl also denotes the summation Ap™ + Ap°"* has no such a complicated dependence 
and is determined only by the temperature T^q 



B. For steady state 



Let us proceed to searching the definition of the SST pressure. Because photons behave as collisionless particles, 
the relation (0) holds for the steady state with replacing the temperatures appropriately for each term (and with 
recalling the closed surface C), Apj_ (T™, To^t) = Ap|(To„f) + ( Ap™(Ti„) + Ap°"*{Tout) ), where the dependence on 
the temperature is explicitly presented. Then via eq. ||2Jl and we obtain 
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2 ( Peqmn) + Peq{Tout) ) ~ ( Apf (T„ ) - Ap^"* (T,,,,) ) 
( T'ln + T^ut ) —Jhk ■ 



(5) 
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Because it is appropriate to define the pressure using the momentum flux even if for the cases of general non-equilibrium 
systems, we define the SST pressure P referring to Apj_. Here, respecting the spirit of the SST mentioned in section 
n we should define P with extracting the effects due to Jhk from Apj^, 



P=— iT'^ ) 

^ — 2^ V «n ^ out ) 



(6) 



Note that this SST pressure reduces to Pf,q in equilibrium cases Tin — Tout, and that P is equal to the average 
of two equilibrium pressures, P — [Peq{Tin) + Peg (Tout ) Further it should be remarked that the value of P 
can be determined operationally via the following steps. First measure the total momentum Ap j_ (Ti„ , Tout ) . Next 
interchange the temperatures of black bodies Tin Tout and measure the total momentum Ap^(ro„t, Tj„). Finally 



the SST pressure is known by calculating their average, P = [Ap^(Ti„,Toi 



Apj_(To„t,r„)]/2. 



For the end of this section, it is necessary to comment that the eq. Q may lead one not to eq. (O but to 



Apj. {Tin, To 



Peq{Tout) + (2(t/3c — Q;) Jfik/uAin- If wc uscd this form of Ap±{Tin, Tout), the SST pressure would 



have been defined as P = Peq{Tout) by extracting the term of Jhk- We have rejected this way of defining P as an 
ill-defined one, since the dependence on the temperature Tin has disappeared. 
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FIG. 3: Two identical steady states. The region B is not measured into the SST volume. 

III. SST FREE ENERGY 

Taking the SST pressure P as the basis, we proceed to searching the definitions of the other state variables 
appropriate to the SST for the radiation field. Further as mentioned in section ^ the extcnsivity/intensivity of state 
variables and the convexity/concavity of the free energy do also play the role as the basis for searching the state 
variables for the steady states. So for the first we discuss about the scale change of the system to make the meaning 
of the extensivity/intensivity in the SST be clear. Then we suggest the definitions of the free energy for the SST (SST 
free energy) and the other state variables. 

A. Scale change and the Extensivity/intensivity 

Someone may think that the scale change is a similarly enlargement of the inner black body and the cavity together 
by the same enlargement rate, or others think it as a similarly enlargement of only the inner black body or the cavity. 
The scale change of the system we consider here includes not only both ways of the similarly enlargement, but also 
arbitrary deformation of the geometrical shape of the inner black body and the cavity. Further, as is the case of 
equilibrium states, all steady states of the radiation field which have the same volume but have different shape of the 
system are identified as the same steady state. Here we should be careful in measuring the volume for the SST (SST 
volume). The SST volume should be measured with the region where both photons emitted at the inner black body 
and at the cavity exist together. For example see figure 01 the region B does not contribute to the SST volume, since 
photons are coUisionless particles. The effect of the scale change on the steady state is simply the change of the SST 
volume V 

Scale change : V XV ( A > ) . 

The extensivity and the intensivity of state variables are defined under this scale change. If a quantity Q is invariant 
under an arbitrary scale change, Q is intensive. If Q changes to \Q, it is extensive. The SST pressure P defined by 
eq. © is obviously the intensive state variable, and the SST volume V is extensive by definition. We assume that all 
the state variables in the SST are classified into to two categories, intensive variables and extensive variables. 

B. Temperature difference and SST temperature 

Before the discussion of the SST free energy, we search the state variable appropriate to specify the non-equilibrium 
order of the steady states of the radiation field. One can recommend two candidates for such a variable, the house- 
keeping heat rate Juk and the temperature difference r = Tin — Tout- However it is obvious that J^k is neither 
extensive nor intensive, since Jhk oc Am- Therefore we reject Jhk I'^ contrast r is obviously intensive. Further 
r = for equilibrium cases, which is the plausible property for specifying the non-equilibrium order as mentioned in 
section n Hence we adopt the temperature difference r as the intensive state variable for specifying non-equilibrium 
order. 
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Next let us suggest the definition of the temperature for the SST (SST temperature) T. We require T to be 
intensive, and to reduce to the ordinary temperature in equilibrium cases. Then we define as 

r = i (r„, + To„t) . (7) 

This temperature never takes negative value. One may think it is abrupt to give this definition here. As is going to be 
shown below, this definition of T together with P, r and V leads to consistent definitions of the other state variables. 



C. SST free energy and the other state variables 



From above, the SST pressure P is expressed as a function of T and r. 



T + 5 



T - - 
2 



Since the SST's spirit requires that the same mathematical formalism as the equilibrium case holds for steady states, 
it is appropriate to assume that the SST volume V is the extensive state variable which accompanies the SST pressure 
P. Then we require for the SST free energy F to satisfy the following relation, 

dF_\ 



PiT,r)^-{ — 



This requirement denotes that F can be treated as a function of T, r and V . It should also be required that F{T, r, V) 
reduces to the ordinary free energy, Feq ~ —{Aa/3c)T*^V, for equilibrium cases. Further as mentioned in section 
n we require that F retains its properties of equilibrium cases about the curvature, that is, F{T, r, V) should be 
convex about extensive variable V, d^F/dV^ > 0, and concave about intensive variables T and r, d^F/dT^ < and 
d^F/dr^ < 0. Hence we define as 



r + 5 



T- - 
2 



V. 



(8) 



Note that F is obviously extensive, and that this definition of F is equal to the average of two equilibrium free energies, 
F ^ [F,,{T,„,V) + Fe,iTout,V)]/2. 

Then we can introduce the definition of the entropy for the SST (SST entropy) S by 



S{T, T, V) 



dF 
df 
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r^3 
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(9) 



This S is obviously extensive. Note that this S reduces to the ordinary entropy, Seq{Teq,V) = {Wa / 3c) T^^ V , for 
equilibrium cases, and that S is equal to the average of two equilibrium entropies, S = [ Seq{Tin, V) + Seq{Tout, V) ]/2. 
Further we define a non-equilibrium order parameter as the extensive variable accompanied by r. 



^{T,T,V) 



dF 

a7 



T.V 



4ct 
3c 



V . 



(10) 



Note that this ^' reduces to zero for equilibrium cases, and that 5* is given by the difference of two equilibrium 
entropies, * = [5eg(r„„ V) - Seq{Touu V) ]/4. 

Let us comment about the chemical potential here. In equilibrium cases, the chemical potential for the radiation 
field is exactly zero, since the number of photons are not an independent variable but determined by the volume and 
temperature. In our steady state, since photons are coUisionless particles, the number of photons supplied by the 
inner black body is not affected whether Tin equals Tout or not, and the same is applied to the photons supplied by 
the outer black body. So we require that the SST chemical potential is also exactly zero, fi = 0. 
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IV. SST GIBBS-DUHEM RELATION 

We proceed to searching the definition of the internal energy for the SST (SST internal energy) E. As is required 
for the other state variables in the previous section, we require that E is extensive and reduces to the ordinary internal 
energy, E^q — (Aa/c) T^^ V, for equilibrium cases. Further remember that, for equilibrium cases of radiation field, E^q 
is related to the equilibrium free energy by E^q = —3Feq. Then we define E as 



2cr 

E = -3F^ — 
c 



V. (11) 



Note that this is equal to the average of two equilibrium internal energies, E = [EeqiTm) + Eeq{Tout)]/'^, and that 
the following relation can be easily found, 

E = T S - PV + . (12) 

Here it is necessary to discuss about the variables which E depends on. In equilibrium thermodynamics, E^q 
depends on the extensive variables, Eeq{Seq, V), as is guaranteed by the Gibbs-Duhem relation. Then we can obtain 
the following relation as the SST Gibbs-Duhem relation using P, S and ^' defined in previous sections, 

-SdT + VdP-'^dT = 0. 

Indeed via eq. H12|l . we find the relation, 

dE = TdS -PdV + Td^ . 

This implies that the SST first law holds, and denotes that E is the function of S, V and ^' to lead the foUowings, 

Finally let us point out that the Legendre transformation holds, 

F{T,T,V)^E{S,-^,V)-TS-T^ {=~PV). 

The above denote that our definitions of E and F are consistent as the steady state thermodynamic functions from 
the view point of the functional relation and the dependence on the state variables. 

A supplemental discussion supporting the validity of the definition Hll(l is placed in the appendix, where the total 
energy carried by all photons are considered and the effect of J^k is extracted from the total energy to give the form 



V. SUMMARY AND DISCUSSION 

Respecting the spirit (hypothesis) of the SST that the same mathematical formalism as the four laws of the 
equilibrium thermodynamics exists for the steady states after extracting the effects of the house-keeping heat from 
various quantities, we have searched consistent definitions of the steady state thermodynamic functions and the other 
state variables for the radiation field sandwiched between black bodies of different temperatures. It should be noted 
that the value of the SST pressure P can be operationally determinable, and that all the other state variables we have 
defined for steady states can also be determined operationally, since their definitions are based on P. The remarkable 
results are that the steady state thermodynamic functions F and E are related by the Legendre transformation, 
and that they retain the plausible properties, the convexity /concavity and the extensivity. These facts give us the 
expectation that the SST does exist at least for the radiation field. The most important and essential issue remained 
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to be settled is the SST second law. In the next work, we should check whether the S defined by eq. @ never 
decreases or not along any change of steady states. 

Next let us discuss about the symmetry under the interchange of the temperatures, t — r. Because the state 
variables we have defined for steady states are related to the average or difference of the state variables for equilibrium 
cases, it is easily found that r and 4* are anti-symmetric and all the other state variables are symmetric. In other 
words, the state variables which vanish in equilibrium cases are anti-symmetric, and the state variables which survive in 
equilibrium cases are symmetric. This symmetry together with the relation to the equilibrium variables (the average 
or the difference) may be a guiding principle to define state variables for the steady states for the radiation field 
produced by many black bodies of different temperatures. For example, the SST pressure for the system composed 
of three black bodies may be P = [Peq{Ti) + Peq{T2) + Peq{T2) ]/3. 

Finally we discuss what the situation shown in figure |2| (right side one) imphes. Let Vsst and Vb denote the SST 
volume and the volume of the region B respectively. Take the equilibrium limit r — > 0. Then the SST entropy 
becomes S (16(t/3c) Tg^ T4st- On the other hand, in this limit, the radiation field in the whole region of volume 
Vsst + Vb is in the equilibrium state of entropy (16cr/3c) T/^ [Vsst + Vb) = S{Teq,T = 0, Vsst) + Seq{Teq, Vb). Hence, 
it may be implied that the simple additivity of the SST entropy and the equilibrium one holds for the system where 
a steady state coexists with an equilibrium state. If this is the case, the total entropy of the system shown in figure^ 
may be given by, S + Sm + Sout, where S is the SST entropy of the radiation field, Sm is the equilibrium entropy of 
the inner black body and Sout is that of the outer one. Further if the SST second law holds with this additivity, the 
relaxation processes to a total equilibrium state of our system composed of two black bodies can be traced within the 
phenomenological framework, which can be applied for studying the black hole evaporation. 
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APPENDIX A: SST INTERNAL ENERGY 

We consider the concentric spherically symmetric system where the radii of the inner black body and the cavity are 
Rin and Rout respectively. Then we derive the definition of the SST internal energy using a mean total energy 
carried by all photons. First define the total mean energy (e) by 

{e)^aT,lA„ m)+aT^i,Aout {ti) , 

where Ain = AttR^^ and Aout = ^''^Rout- The quantity is the number average of flight time of photons emitted at 
a point r on the inner black body's surface, and {t[) is that of photons emitted at q on the cavity's surface. 

Let a be the angle between the direction of a photon's orbit and the normal direction to the inner black body's 
surface, then the number of photons n{uj) emitted at r par a unit area, a unit time, a unit frequency and a unit solid 
angle is 

-i^exp(/?rM-l'°'"' 

where /3i„ = l/ksTm and uj is the photon's frequency. The flight time ti^{a) of a photon emitted in the direction of a 
is given by 

cit(a) = J - sin^ a - i?i„ cosa . 



FIG. 4: Concentric spherical system. 



Therefore the averaged time is given by 



— — J dsdujn{uj,a)ti^{a) 

9 / P \ 9R / p2 \ 3/2 



3Ri 3 3c \ R. 
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where ds is the sohd angle element and N — j dsdLun{uj,a). Next let b be the angle between the direction of a 
photon's orbit and the normal direction to the cavity's surface, then the flight time t|(a) of a photon emitted at q in 
the direction of b is given by 



'2Rout cos 6, for b^ <b < 7r/2 

Rout cos b - \ RL - Ro„f sin^ b, for < 6 < 5n 



ctiib) = 

■ 'tout COS - Y it -„ - n-^f. sm- o , lor u <. o <_ oq 
where 6o is shown in figure 21 Therefore the averaged time is given by 



C \ 3 3-Rq„( y 3c \ -Rout 



2 \ 3/2 



Consequently we obtain 



3/2 



where V = (47r/3) (i?^^^ - R^^). Then, according to the SST's spirit, we define the SST internal energy E by 
extracting the term of Jhk, 



This is the same form as eq. 
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